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A non-reciprocal phonon model for microwave or optical isolators is considered. It gives a simpler 
framework to further investigate the formerly argued possibility for a heat transfer between black 
bodies at common temperatures. While the non-dissipative device retains the Detailed Balance 
property, the presence of dissipation breaks it. This property allows a net transfer of heat between 
the two black bodies at common temperatures, whenever the absorptive elements are at lower 
temperatures than the one being common to the bodies. 
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In this work we continue the previously started exam- 
ination of the physics of the heat transfer among bodies 
and the fact that such energy flows are expected to al- 
ways spontaneously occur from higher to lower temper- 
ature regions. This property is equivalent to the Law 
of increasing Entropyu. However, in our view its theo- 
retical derivation from statistical physics deserves more 
close attention. A concrete fact underlining this need 
is the lack of definition of the entropy for physicaLpro- 
cesses occurring far from the equilibrium situationo. It 
can be understood that for a very wide class of physical 
systems having reciprocal dynamics, the heat will always 
travel from higher to lower temperature regions. How- 
ever, when the bodies have a non-reciprocal dynamical 
response, this property is no so naturally evident. Non- 
reciprocal phenomena occurs in practice witHa-deiiices 
such as optical and microwave isolators andotffll3'Eland 
more generally in all the physics of the magnetism hav- 
ing so many applications in technology. 

The former paperi made an analysis of two types of 
isolators connecting two regions containing black-body 
radiation at common temperatures. In the present article 
we consider the realization of the same kind of phenom- 
ena within a simple soluble ID non-reciprocal phonon 
model. It is defined by a string having a uniform charge 
density which is placed within a magnetic field orientated 
along the equilibrium line of the string. The systems 
wave modes are explicitly determined. They show two 
independent Faraday waves for each direction of propa- 
gation. These are employed for defining a phonon version 
of the optical and microwave isolators discussed inS. 

A main new possibility introduced here is the addition 
of a lateral string (the "vane") at one of the ends of the 
isolator which simulates the dissipative resistance vane 
in the microwave or optical versions of these devices. It 
is then shown that when the vane is disconnected, the 
transmission coefficients through both sides is identical 
no matter the frequency or the orientations of the slits 
simulating the optical polarizers. This result in fact rep- 
resent a clear .solution of the so called Wien paradox in 
optical systemB. Also it illustrates the validity of the De- 
tailed Balance principle for the two terminal conserva- 



tive version of the isolator. After that, the waves modes 
for the system including the vane were also solved. The 
calculated spectrum for the transmission coefficients im- 
plied the breaking of the Detailed Balance principle for 
the isolator including the vane. However, it follows that 
when the three outputs are at the same temperature, the 
net flux through any of the three terminals exactly van- 
ish. This result excludes the presence of the heat flux 
unbalance under full thermal equilibrium. However, a 
net heat flow between Black-Bodies at common temper- 
atures is predicted when the vane is placed in contact 
with a very low temperature heat reservoir. 

Seeking for a simple model -for the optical and mi- 
crowave isolators discussed inlj let us consider a La- 
grangian system formed by a uniformly charged string 
in presence of an also uniform magnetic field. The La- 
grangian of the system and its Euler equations are given 
by 

L = J dx{^{d t Hf-^(d x -ltf + 

(1) 

= - Pm d 2 lt + k h d 2 x lt + ^-{d t d x lt) x ~S. (2) 

where lt(x, t) are the vector displacement of the points 
of the string with respect to their equilibrium positions 
and p m and p q are the linear mass and charge densities, 
kh is the Hook law constant. The magnetic field is ori- 
ented along the string in the x direction B — (1,0,0). 
It should be mentioned that the self interaction of the 
charged string will be disregarded in order to simplify 
the discussion. The coulomb interaction, basically will 
produce the appearance of a plasmon mode introducing 
a gap at low frequencies. Searching for the oscillation 
modes in the usual form 

Ht(x, t) = exp(i.kx — ik x )lt(k) 

in which k Q = w/c s ,x — c s t and c s = \JkhJPm is the 
sound velocity in the cord, c is the light velocity. The 
Euler equation in (1) after written in components takes 
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the form 

( (fc2 - ^ + lI ^^ B ) * = °' 

which after to be solved produce the dispersion relations 

for longitudinal and transversal waves for positive energy 

values e = hk D c s k a = k,k a — oM + \Jk 2 + A/ 2 , a = 
± M= j*£|c. 

Inversely, the momenta and frequencies obey the rela- 
tions k a — (k a — oM) 2 — M 2 , with a = ±. It should 
be noticed that for a fixed energy e = h k a c s there are two 
solutions having momenta k + and fc_ . After perform- 
ing linear combinations of two modes associated with the 
same fixed positive frequency, Faraday waves can be ob- 
tained. That is, they show the Faraday effect: their po- 
larizations are linear at any point and its angular orienta- 
tion changes proportionally with the displacement along 
the string. 

A phonon version of the experimentally existing mi- 
crowave or optical isolators investigated irjj is accom- 
plished by considering an arrangement depicted in Fig. 
1. 




FIG. 1: Phonon Isolator 

The device can be described in the following way: a) 
Two slits forming an angle 9 among them are situated 
a distance L restricting the movement of a long piece 
of the string described in previous section. They lay 
within a plane being orthogonal to the string and it is 
assumed that there is no friction between them and the 
cord. In other words they form conservative mechanical 
constraints on the system, b) At a distance S <C L of 
one of the them, let us say the right one, a second string 
is attached to the main one, by mean of a bar sliding 
without friction inside a cylinder at rest. That is, the 
point at the junction is constrained to move linearly in 
a direction orthogonal to the initial string. This fact im- 
plies that the newly added string, that we will call from 
now on, the "vane", is allowed to oscillate only longitu- 
dinally, c) Let us also consider that the frequency lies in 
the linear part of the dispersion relations, that is, when 



M <C k Q . Then, the length L is chosen to assure that a 
Faraday mode propagating to the right and starting at 
the slit 1 being polarized along it, after arriving to the 
slit 2 will be also polarized along it. This property, in 
turns will assure that the waves incident on the left of 
slit 1, being polarized along it will have a transmission 
coefficient equal to one, whenever the coupling with the 
vane is disconnected, d) Next, the angular orientation of 
the vane at the fixing point distant a length 5 from the 
slit 2 is selected to coincide with the exact polarization 
of the before analyzed Faraday wave having transmission 
coefficient equal to one. e) At last, let us consider that 
outside the two slits, the charge density of the string van- 
ish and also that additional non-dissipativc mechanical 
constraints exist which only allow the transversal polar- 
ization of the waves in the directions of the associated 
slits. The only purpose of this assumption is to reduce 
the number of independent wave modes of the system in 
order to simplify the discussion. It should be notice that 
in the proposed arrangement the longitudinal oscillations 
do not couple with the transversal ones (within the linear 
approximation only). 

The above conditions, then assure that the system can 
be described by a three terminal scattering matrixEj: 
h = Sijdj, i,j — 1,2,3, where the pair (ai,fei) con- 
tains the incident a± and reflected b\ wave amplitudes 
at the slit one and (0,2, 62), (oi, &i) are the corresponding 
quantities at the slit 2 and the vane respectively. 
In order to solve for the matrix S, auxiliary Faraday 
modes propagating in the two internal regions, shown in 
Fig 1 were defined. Then, the amplitudes for such modes 
were related among themselves by imposing the continu- 
ity of the amplitudes and the components of the forces 
along the directions of the slits 1 and 2 and the vane's 
direction. Further, after defining the scattering modes 
associated to the input, let say the one with number 
i(i = 1, 2, 3), as such wave modes for which the incident 
amplitude a,j = 1 and all the other incident amplitudes 
Oj =0, j ^ i at the rest of the terminals vanish, the set 
of linear equations defined by the boundary conditions 
were solved. More precisely, as one of the propagating 
modes shows a gap for low frequencies, this wave turns 
to be spatially damped boundary mode when the fre- 
quency is low enough to be in the gap. The effects of 
these low frequency oscillations were disregarded in this 
initial discussion. 

The above described equations for finding the scatter- 
ing modes were solved initially for the case in which the 
vane is disconnected from the central string. This situ- 
ation corresponds to an optical or microwave isolator in 
which the " resistive element" absorbing the power in the 
output 2 is absentQ. The system is also completely analog 
to the one related with the so called Wien paradoxel. The 
transmission coefficient for the resulting two terminal de- 
vice T = j^j, 0,2 = 0, for the wave incident upon the slit 
one was then evaluated as a function of the wave vector 
and the angle between the two slits. The results are 
shown in the two dimensional plot in Fig. 2. It, should be 
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reminded that the transmission coefficient for the wave 
coming in through the slit 2 can be extracted from the 
figure by changing the sign of the angle 9. The data cor- 
respond to a length L for which the angle 9 — tt/4 pro- 
duces a value one for the transmission coefficient from 
left to right at high frequencies. The picture remarkably 




FIG. 2: Transmission coefficients for the Phonon Isolator 
without losses 

shows that in the absence of the vane, in spite of the non- 
reciprocity of the system, the transmission coefficient as- 
sociated to the two inputs are exactly equal whatever the 
value of the angle and the frequency. In other words, the 
fraction of the incident power transmitted from one side 
to another is exactly the same whatever side is taken as 
the input. The equality to one of the transmission for the 
angle 8 = ? is also evidenced. The results curiously indi- 
cate that in spite of the multiple reflections between the 
two slits of the wave entering through slit 2, the trans- 
mission coefficient for this wave turns to be equal to one. 
This conclusion gives-,a clear solution of the above men- 
tioned Wien paradoxp and also implies the exact validity 
of the Detailed Balance principle for the non-dissipative 
isolator. In conclusion, the absence of a net heat flux 
between two black-bodies joined by such isolators also 
follows. 

Further, the case in which the vane is connected was 
investigated. At this point it may be useful to remark 
that we considered in detail the region of frequencies 
M <C k Q in which the rotation of the polarization per 
unit length is a constant and its value for the length 
L was selected to produce the transmission coefficient 
T = 1. Therefore, for this situation, we evaluated the 

\b 2 

two transmission coefficient £12(^0) = a * 2 ; a i — a 3 — 0: 

16 2 

and ti3(k ) = -j^p- , al = a 2 = 0, in the now three termi- 
nal device and for the angle 9 between the slits equal to 

7T 

4 ■ 

As the reflection coefficient for the wave incident at the 
slit 1 is equal to zero, the question about a possible power 
flow unbalance at the input number one is determined by 



these quantities. The values of the transmission coeffi- 
cients are illustrated at Fig. 3 as functions of the energy 
in the range k /M 6 (3, 10) where the linearity condition 
of the dispersion relations is obeyed very approximately. 

Figure 3 clearly indicates a breaking of the detailed 
Balance in the three terminal system simulating the op- 
tical and microwave isolators. This conclusion follows 
from the fact that £12(^0) 7^ £21(^0) = 1- The possibility 
for the existence of physical systems not satisfying this 
property was elderly advanced by the same BoltzmannEI 
long time ago. 

However, let us now argue that under the assump- 
tion of the thermal contact of the three terminals with 
heat reservoirs characterized by a common value of the 
temperature, the net heat flux through any of the in- 
puts is equal to zero. For the slit number one this out- 
come is a direct consequence of the results in Fig. 3 
and the identical formula for the link between the flux 
of energy and the amplitudes of the waves in each of 
the outputs. This last property is determined by the 
assumed common physical parameters (density p m and 
Hook law constant kh ) of the strings associated to any 
of the terminals. The relevant property of Fig. 2 de- 
termining the vanishing of the heat flux at the left of 
the slit 1 is the relation ti 2 (k ) + £13(^0) = 1, which is 
clearly satisfied by the data in the picture. This is also 
the situation for any of the other terminals as a general 
consequence of the consepative character of the three 
terminal device (See Ref.EEl), which implies that the S 
matrix is unitary. To see this explicitly, it is useful to 

consider that the columns of S defined by a fixed value of 

(7) 

j, give the components b\ f , i — 1,2,3, of the scattering 
modes normalized by the condition ctj = 1 and a/ = 
for I ^ j. Therefore, let us consider the power 
the waves emerging from output % under the assumption 
that equal powers are entering through the three termi- 
nals due to the common temperature conditions. Then 
it follows P e (,;) = Po £j=i bf(hf)* = p £j=i Sq S% = 
Podia* = p where p is the constant relating the square 
of the amplitudes (defined by some convention) with the 
associated power flow in a travelling wave. The above re- 
lation implies that the power of the outgoing waves at any 
of the terminals is identical to the power of the thermal 
incident waves entering through the same input, when- 
ever all the heat reservoirs are at a common temperature 
T. 

Finally, we will consider that the vane terminal is lay- 
ing at a very much small temperature T v « T where 
T is the temperature of the left and right heat reservoirs 
in contact with the pieces of the central string arriving 
at the left and right slits 1 and 2. Also the angle 9 will 
be selected to furnish a transmission coefficient equal to 
one for the waves incident through the slit 1. Under this 
conditions it follows that at least for a very extended re- 
gion of the energies for which the relation M <C k a is 
valid, the net flux passing from left to right of the iso- 
lator is greater than the one circulating in the opposite 
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FIG. 3: Transmission coefficients determining the outgoing 
heat flow in the input at the left of the slit 1 



sense. Since this property is independent of the temper- 
ature, we will assume that it can be elevated sufficiently 
enough for to be able for disregarding the low energy re- 
gion, when deciding about the question of the net heat 
balance. Under this central assumption, it can be con- 
cluded from the data in Fig. 2 that the system will pro- 
duce a heat flow passing from the Black -Body at the 
left of the slit 1 to the black body at the right of slit 2. 
This conclusions motivate to consider in further works 
the possibility for a reduction of the statistical entropy 
through concrete mechanisms for the extraction of energy 
through the vane. In ending, I would like to acknowledge 
the discussions with Dr. A. Gonzalez, C. Trallero and F. 
Comas and the general support of Dr. D. Sheehan, the 
Org. Committee of the Conference and the Associate 
Scheme of the ASICTP. 
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